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ABSTRACT 

Spectral synthesis in 3-dimensional (3D) space for the earliest spectra of Type la supernovae 
(SNe la) is presented. In particular, the high velocity absorption features that are commonly seen 
at the earliest epochs 10 days before maximum light) are investigated by means of a 3D Monte 
Carlo spectral synthesis code. The increasing number of early spectra available allows statistical 
study of the geometry of the ejecta. The observed diversity in strength of the high velocity 
features (HVFs) can be explained in terms of a "covering factor" , which represents the fraction 
of the projected photosphere that is concealed by high velocity material. Various geometrical 
models involving high velocity material with a clumpy structure or a thick torus can naturally 
account for the observed statistics of HVFs. HVFs may be formed by a combination of density 
and abundance enhancements. Such enhancements may be produced in the explosion itself or 
may be the result of interaction with circumstellar material or an accretion disk. Models with 1 
or 2 blobs, as well as a thin torus or disk-like enhancement are unlikely as a standard situation. 

Subject headings: supernovae: general — radiative transfer — line: profiles 



1. Introduction 

Type la supernovae (SNe la) have found use in 
cosmology after it has been established that the 
maximum luminosity can be calibrated via an em- 
pirical relation between it and the shape of the 
light curve (LC, Phillips 1993; Riess, Press & Kir- 
shner 1996). The origin of this empirical relation 
is however not fully understood (e.g., Mazzali et 
al. 2001), mainly because of uncertainties regard- 
ing the properties of the explosion mechanism. 

A one-parameter ordering scheme similar to the 
one found for the LC was also proposed for the 
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spectra (Nugent et al. 1995). In this scheme tem- 
perature is the driving parameter characterizing 
the spectra. While this relation may apply near 
maximum, detailed observations of SNe la with ex- 
tended early time spectral coverage indicate that 
SNe la with similar LCs may have different absorp- 
tion line velocities. In particular, the absorptions 
of Si II A 6355, SUA 5640, Ca H H & K and Ca 
II IR triplet have been investigated by several au- 
thors (e.g., Patat et al. 1996; Hatano et al. 2000; 
Kotak et al. 2005; Benetti et al. 2005). 

Simultaneously, the number of polarization ob- 
servations has increased (SN 1996X, Wang et al. 
1997; SN 1999by, Howell et al. 2001; SN 2001el, 
Wang et al. 2003; SN 2004dt, Wang et al. 2004) 
and it has become possible to infer asphericity by 
means of spectropolarimetry. These results sug- 
gest that some SNe la have a somewhat aspher- 
ical geometry of the photosphere (> 10% for SN 
2001el or ~ 20% for SN 1999by, assuming a ellip- 
soidal geometry). Even larger degree of aspheric- 
ity is inferred for the distribution of intermediate 
mass elements (at the level of < 25 %; see Wang 
et al. 2003). Multi-dimensional numerical simu- 
lations have also suggested that the explosion is 



asphcrical, mainly owing to the nature of the de- 
flagration flame in 3D numerical simulations (Rei- 
necke, Hillebrandt fc Niemeyer 2002; Gamezo et 
al. 2003; Ropke & Hillebrandt 2005). 

Recently, high velocity features (HVFs) in the 
Ca II IR triplet have been the subject of inter- 
est: SN 1994D (Hatano et al. 1999), SN 1999ee 
(Mazzah et al. 2005a), SN 2000cx (Thomas et 
al. 2004), SN 2001el (Wang et al. 2003), SN 
2003du (Gcrardy et al. 2004), SN2005cg (Quimby 
et al. 2006). Here, HVFs are defined as absorp- 
tions with velocities much higher than the pho- 
tospheric component, which has a typical veloc- 
ity V ~ 15,000 km s~^ at ~ 1 week before max- 
imum. HVFs are often seen detached from the 
photospheric component, at velocities of 17,000- 
29,000 km s~^. The above papers discuss the pos- 
sible origin of HVFs including primordial metallic- 
ity (Lcntz ct al. 2000), a property of the explosion 
and interaction of the SN ejecta with circumstellar 
matter (GSM). Whatever their origin, understand- 
ing HVFs can cast light on the explosion mecha- 
nism. 

Mazzah et al. (2005a) tried to fit the HVFs of 
SN 1999ee using abundance or density enhance- 
ments. Their results suggest that it is impossi- 
ble to reproduce HVFs by abundance enhance- 
ment only and ~ 0.1 Mq of additional material 
is required if only the density enhancement is as- 
sumed. Gerardy ct al. (2004) and Quimby et al. 
(2006) discussed how HVFs may result from in- 
teraction with a GSM. They obtained good agree- 
ment for several features including Mg II, and 
found 2 X lO^^Afg and 5 - 7 x IO'^Mq of the 
GSM are needed for SN 2003du and SN 2005cg, 
respectively. It should be noted that the explo- 
sion model (i.e., abundance distribution and the 
density structure) and the way to introduce addi- 
tional material are not identical with each other. 
Recently, Mazzali et al. (2005b) have shown that 
almost all SNe la have Ga II IR triplet HVFs at 
the earliest epochs. These HVFs have velocities 
ranging from 17,000 km s^^ to 29,000 km s~^ in 
different supernovae. The ubiquity of HVFs may 
indicate that they do not come from an extreme 
environment. Probably, the combination of differ- 
ent phenomena results in the observed HVFs. 

The geometry of the ejecta has been studied in 
multi-dimensional space. Kasen et al. (2003) an- 
alyzed both spectroscopy and spectropolarimetry 



of SN 2001cl and showed that both an asphcrical 
photosphere and a single high velocity blob can 
reproduce the observations. Thomas et al. (2004) 
derived a similar geometry for the Ga II HVFs of 
SN 2000CX. 

Although all previous studies performed mod- 
eling for each SN, no systematic study has been 
made of the 3D properties of the geometry. We 
present here synthetic spectra computed in 3D 
space for various gcxnricitries and show how HVFs 
are affected by different geometrical configurations 
and line-of-sight effects. The increasing number of 
early spectra available enables us to investigate the 
statistical properties of HVFs, and to constrain 
the geometry of the ejecta and possibly the na- 
ture of the explosion. 

In §2 we discuss how HVFs are formed and sug- 
gest a possible parametrized description. In §3 we 
present the method of calculation for the spectra 
and our models. In §4 the results and some prop- 
erties of 3D computations are presented. In §5 
the probable geometries that can reproduce the 
observed trend are discussed. Finally, conclusions 
including considerations on the origin of the high 
velocity material and on the asphericity of the ex- 
plosion are made in §6. 

2. Properties of the HVFs 

HVFs come in different forms. Figure 1 shows 
the Gall IR triplet in the earliest spectra of SN 
1999ee (-9 days before B maximum), SN 2002dj 
(-11 days), SN 2003cg (-8 days), and SN 2003du 
(—11 days). At that epoch, the photospheric ve- 
locity is Vph ~ 12, 000 - 14, 000km s"!. The pho- 
tospheric component of the Gall IR triplet has 
V - 14,000 - 16,000km s^^ in all SNe la, as 
suggested by the data of Mazzali et al. (2005b). 
This is slightly higher than the photospheric veloc- 
ity, since the Gall IR triplet is stronger than Sill 
or SII, and forms above the photosphere. How- 
ever, the velocities of the absorption minima in SN 
1999ee (- 20,000 km s'^), SN 2002dj (- 25,000 
km s-i) and SN 2003du (~ 20,000 km s'^) are 
much higher, suggesting that these lines form well 
above the photosphere. Therefore the HVFs are 
expected to have a different origin that the photo- 
spheric absorption. Gomparing SN 2002dj and SN 
2003du, which were observed at the same epoch, 
it can be noticed that their HVFs have different 
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strength and velocity range. The same is true for 
SN 1999ee and SN 2003du, although their spectra 
were not taken at exactly the same epoch. At very 
early epochs, even a small difference in epoch can 
strongly affect the strength of the photospheric 
component. However, judging from a number of 
early spectra we can safely conclude that there are 
real variations in the strength of the HVFs. There- 
fore we examine a parameterization of the prop- 
erties of the HVFs based on their velocity range 
{vuv), their strength (-S'hv) and on the strength 
of the photospheric component (Sph). 

We performed numerical tests to verify what 
governs the properties of the HVFs. The Monte 
Carlo SN spectrum synthesis code described in 
Mazzali & Lucy (1993), Lucy (1999) and Maz- 
zali (2000) was used. The code computes a syn- 
thetic spectrum based on the luminosity, the po- 
sition of the photosphere, the time since the ex- 
plosion, the density structure and the abundance 
distribution. The density structure of the stan- 
dard deflagration model W7 (Nomoto, Thieleman 
& Yokoi 1984) was used in all the calculations 
presented in this section. Photospheric velocity 
{Vph = 12,500km s~-^), epoch (8 days since the 
explosion assuming a risetime of 19 days) and lu- 
minosity (L = 3.0 X 10^^ erg s~^) were fixed for 
simplicity. 

First we consider the case where Sph and whv 
are fixed and parameterize S'hv through the line 
optical depth of the CaH IR triplet. We as- 
sume that the high velocity region lies at 77hv ~ 
22,000 - 29,000 km s^^. In this region, the line 
optical depth is imposed irrespective of the value 
compiited consistently with the W7 structure in 
order to simplify the investigation of the behav- 
ior of HVFs. For example, if line optical depth is 
taken as 1 — exp(— r) = 0.5, about half the photons 
that come into resonance with the Call IR triplet 
are scattered. Figure 2 shows HVFs with various 
line optical depths 1 — exp(— r) = 0.1 — 0.9. With 
this approximation the depth of the HVFs depends 
almost linearly on 1 — exp(— r). The origin of the 
observed diversity in strength is investigated in 3- 
dimensional space in §3 and §4. 

Next we show how the strength of the photo- 
spheric component affects the absorption feature. 
Here Sph is parameterized through the abundance 
of Ca. In the calculation, a homogeneous abun- 
dance distribution is used, which is a reasonable 



assumption because the region used in the com- 
putation is only v > 12,500km s~^. Figure 3 
shows the Call IR triplet with various strengths 
of the photospheric component. The mass frac- 
tion of Ca is X(Ca) = 0.068,0.034,0.017,0.008, 
and 0.04, respectively. The photospheric absorp- 
tion depth does not change linear with the abun- 
dance, and the line becomes saturated if the abun- 
dance is more than ~ 1%. Since in this calculation 
the line optical depth at high velocity is fixed as 
1 — exp(— r) = 0.7 and the velocity separation is 
large, the HVFs are not affected by the photo- 
spheric component. 

Using only the abundance of Ca to change the 
strength of the photospheric component is not 
completely realistic because the line strength also 
depends strongly on temperature and ionization 
state. Ca is mostly Calll at the typical temper- 
ature of SNe la before maximum (T ~ 12,000K) 
and only a small fraction is Call and CalV. The 
strength of the Call IR triplet is proportional to 
temperature because of the behavior of the ion- 
ization. We tested models with different lumi- 
nosities, which yield different temperature struc- 
tures. We found that even with the highest Ca 
abundance (X(Ca) = 0.068), when the luminos- 
ity is increased by a factor of ~ 4, the strength 
of the Call IR triplet becomes weaker than in 
the original model with the lowest Ca abundance, 
X{C&) = 0.002. An inverse correlation between 
the strength of the photospheric component and 
the luminosity at these early epochs might there- 
fore be expected. To veriiy this, more observations 
are needed and it is necessary to fit these individ- 
ually to determine the exact photospheric velocity 
and temperature. This is beyond the scope of this 
work, and will be the subject of a separate inves- 
tigation. (M. Tanaka et al., in preparation). 

The results of the tests we performed by chang- 
ing 5hv and Sph (Figures 2 and 3) show clearly 
that changing these two parameters is not enough 
to reproduce all the observed HVFs. As seen be- 
tween SN 2002dj and SN 2003du in Figure 1 (and 
see Figure 4 in Quimby et al. 2006), the width of 
the high velocity absorption is not identical among 
SNe that have a similar depth of the HVFs. No 
combination of S'hv and Sph, however, can ac- 
count for such a difference. Therefore we conclude 
that the additional parameter vbv is also needed. 
One might think that this is because of our as- 
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sumption of spherical symmetry. However, we ob- 
tain the same results in 3D geometry (§4). Fig- 
ure 4 shows how the three parameters (5hv,i'hv 
and Sph) work. Their combination determines the 
depth and shape of the absorption and the posi- 
tion of the absorption minimum. Although the 
luminosity and the photospheric velocity arc not 
fine-tuned, the shape of the absorption feature is 
reproduced well. 

The numerical tests presented here may be 
somewhat unrealistic because we assumed a de- 
tached spherical shell that is optically thick to 
the the Call line and it is probably more realis- 
tic to consider 3D structures. The explosion may 
be aspherical and has fragmentations (yielding a 
clumpy structure), or interaction between the SN 
ejecta and a circumstellar disk may occur (yielding 
a torus-like structure). In the next section, we use 
3D models to sec what the three parameters used 
in this section (i.e., S'hv,^^hv, and S'ph) physically 
mean. 

3. Method and Models in 3D space 
3.1. 3D Monte Carlo Code 

To investigate 3D effects on the HVFs, syn- 
thetic spectra are calculated assuming various 
geometries (§3.2). We have developed a 3- 
dimensional Monte Carlo radiative transfer code 
based on the ID Monte Carlo code described in 
MazzaU & Lucy (1993), Lucy (1999), and Mazzah 
(2000). The code assumes the Sobolev approxi- 
mation and a sharply defined photosphere above 
which there is no energy deposition. Spherical co- 
ordinates are used in the code in order to treat the 
photosphere precisely. The number of meshes de- 
pends on the model and is typically 20 x 20 x 20 for 
r, cos{9) and 0, respectively. The path of the en- 
ergy packets is traced in 3D space, and the effect of 
back scattering into the photosphere is included. 
Line treatment is free from the approximation 
of resonance scattering, i.e., photon branching is 
considered correctly as outlined by Lucy (1999). 
The free parameters are the luminosity L(9,(j)), 
the photospheric velocity Vp\i{9,(j)) and the epoch 
since the explosion t, all of which are defined as 
in the ID code. Though angle dependent L and 
Vp\i can be used in the code, here we assume them 
to be uniform for simplicity. We again assume 
Vph = 12,500km s~^, t = 8.0 days since the ex- 



plosion and L — 3.0 x 10 erg s~'^as typical values 
for the earliest spectra of SNe la. With these 
parameters, the temperature structure, excitation 
and ionization are computed in all zones using a 
modified nebular approximation (Mazzali & Lucy 
1993). This assumes that there is no net exchange 
of energy between matter and radiation. Photon 
fiux is collected by recording energy packets at the 
outer boundaries of the ejecta and is binned into 
a 50 X 50 solid angle mesh. Since the path of each 
packet is affected by the density and temperature 
structure, the emergent spectrum depends on the 
orientation. 

3.2. Models 

Our models are based on the spherical defla- 
gration model W7 (Nomoto ct al. 1984). We 
introduced additional material in the outer lay- 
ers of the ejecta to produce HVFs. We are not 
concerned here with the origin of this material, 
which may come from fluctuations of the explo- 
sion or from interaction with CSM or an accretion 
disk. Abundance enhancements are not consid- 
ered because they have been shown not to be suit- 
able for reproducing HVFs (Mazzali et al. 2005a). 
The velocity range and the degree of density en- 
hancement are determined by fitting one of the 
strongest Call HVF, that of SN 2002dj, with the 
ID code. As a result, ^ O.IM0 of material at 
V ~ 22, 000 - 29, OOOkm s'^ is added in a sheU. 

We fixed this degree of density enhancement in 
our 3D computations. We used a density enhance- 
ment rather than an artificial line strength, which 
was used in the previous section, because we want 
to explore the cause of HVFs. In fact, a spheri- 
cal density enhancement has the same effect on the 
Call IR triplet profile as the artificial line strength. 
However, other lines, such as SiH A6355, arc also 
affected when using a density enhancement. If the 
observed behavior of both the Call IR triplet and 
Sill A6355 is reproduced consistently with a den- 
sity enhancement, we may conclude that HVFs are 
due to density enhancements. Additionally, the 
mass invoked in HVFs should be reduced in 3D 
space. In Mazzali et al. (2005a), a large amount 
of material (~ O.IOM© in the case of density en- 
hancement only) was required to produce HVFs, 
and this was a serious concern in the hypothesis 
that HVFs come from either the explosion or CSM 
interaction. Although Gerardy et al. (2004) and 
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Quimby ct al. (2006) found only - 2 x IQ-^M© 
and ~ 5 - 7 X 10~^Mq are enough for the HVF of 
SN 2003du and SN 2005cg, respectively, the HVF 
that has a higher velocity and a deeper absorption 
has been observed like SN 2002dj (see Fig. 1). We 
will give an estimate of the amount of material 
needed in a 3D calculation. 

We mapped the density enhancement into 3D 
space conserving the radial velocity range and the 
degree of enhancement derived in the ID calcula- 
tions. We tested various morphologies that may 
be realistic, including one or two large blobs, a 
small number of discrete blobs, a crowded clumpy 
structure, and tori of various opening angles. 

4. Results 

First we describe the general behavior of HVFs 

in 3D models taking a model with a large blob 
(Model Bl) as an example. Figure 5 shows the ge- 
ometry (left) and the synthetic line profile (right). 
In the left panel, the red sphere shows the photo- 
sphere and the blue region shows the density en- 
hancement. The depth of HVFs in the synthetic 
spectra is different for different lincs-of-sight. Seen 
on the z-axis, which is defined as the direction of 
the blob, the absorption at high velocity is deep- 
est. If we move toward the equatorial plane, the 
high velocity absorption becomes shallower and 
shallower and it disappears when the line-of-sight 
reaches the edge of the blob. Note that having 
the blob on one side has only a small effect on 
the emission profile, because the volume occupied 
by the blob is small. The range of profiles in the 
right panel of Figure 5 is qualitatively very similar 
to Figure 2, which shows that our ID parametriza- 
tion captures the role of the blobs. If the blob is op- 
tically thick, the depth of the HVFs is determined 
by the fraction of the projected photosphere that is 
concealed by the dense blob for each line-of-sight. 
We define as the "covering factor" and denote it 
as /. 

Figure 6 shows that the HVFs obtained using 
different covering factors in 3D models resemble 
those obtained with ID synthetic spectra for dif- 
ferent values of 1 — exp(— t). The values are not 
exactly the same because our 3D blobs are not in- 
finitely optically thick. In fact, given a typical op- 
tical depth To in the blobs, one expects a relation 
exp(— r) ~ (1 — /)-|-/exp(— To), where r is the line 



optical depth introduced in the ID computation. 
Thus / = 1 — exp(— r) if tq = oo. Since tq has 
a finite value, / ~ (1 — exp(— r))/(l — exp(— ro)), 
which is smaller than the value in an infinitely 
thick case. Figure 6 clearly indicates this. Com- 
paring the value of 1 — exp(— r) and / that yield 
similar spectra shows that the effective tq of the 
3D blob decreases with increasing angle. This is 
probably because the average length of a segment 
crossing the blob is smaller for larger angles, as 
there are more grazing trajectories. 

If such a blob exists, the variation of the 
strength of the HVFs (S'hv) can be reproduced 
along various lines-of-sight without introducing 
unknown parameters such as line optical depth 
1 — exp(— t) in §2. We roughly classify the HVFs 
by their strength into three groups, strong HVFs 
(/ > 0.7), medium HVFs (/ ~ 0.45 - 0.7) and 
weak HVFs (/ < 0.45). For example, SN 2002dj 
and SN 2003duTiave "strong" HVFs, SN 1999ee 
has a "medium" HVF, and SN 2003cg is "weak" 
(see Fig 1). 

However, we have to consider the statistical 
properties of HVFs. Recent observations suggest 
that a considerable fraction of the earliest spec- 
tra have high velocity absorption. On the con- 
trary. Model Bl suggests that strong HVF will be 
observed only ~ 0.5% in all SNe. Figure 7 is a 
histogram of the covering factor /. It clearly indi- 
cates that only very few orientations can produce 
strong HVFs and about 90% of the observations 
(/ < 0.45) should show only the photospheric ab- 
sorption. Figure 8 shows the observed frequency 
of the HVF's strength made from Mazzali et al. 
(2005a and 2005b) by distributing 8 SNe into 4 
bins. The observed depth of the HVF is converted 
to the covering factor thorough the numerical fit 
as in Figure 4. It should be noted that the num- 
ber of the sample is not enough in Figure 8 but the 
ubiquity of the HVFs is not overthrown even if we 
consider other SNe (Mazzah et al. 2005b). The 
clear differences between the histogram of Model 
Bl (Fig. 7) and the observation (Fig. 8) are seen 
at / 0.0 (too much in Model Bl) and > 0.5 (too 
small in Model Bl). Although it is of course not 
necessarily required that all observations are re- 
produced with a single geometry, a general geom- 
etry may be inferred if a certain model can repro- 
duce the observed trend. We investigate several 
geometrical models in the following subsections 
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and compare the distribution of synthetic prop- 
erties to the observed trend. 

4.1. Blob Models 

We showed that a model with one blob is un- 
likely as a general geometry of SNe la. We there- 
fore vary both the size of the blobs and their num- 
ber. Wc show the expected statistical properties of 
each model by means of the distribution of cover- 
ing factors, since we have shown that there is a di- 
rect correspondence between / and HVF strength. 
Figure 9 shows the distribution of covering fac- 
tors in a model with two blobs (Model B2). As 
expected, the fraction of lines of sight with large 
covering factor doubles and the fraction of small 
factors decreases. This is however still not consis- 
tent with the fact that most spectra have high ve- 
locity absorption. Although we tested various blob 
sizes, most lines-of-sight have / < 0.4 because of 
the large space between the two blobs. We then 
increase the number of blobs adjusting their size 
so that the photosphere is not completely covered. 

Next we consider a case with il blobs. Fig- 

ure 10 shows the distribution of covering factors 
for a model with 6 blobs with opening angle 60 
degrees (top panel; Model B6), a model with 5 
blobs of 80 degrees (middle panel; Model B5) and 
a model with several blobs of 30 degrees (bot- 
tom panel; Model B30). From top to bottom, the 
space between the blobs becomes smaller. Syn- 
thetic spectra show that the correlation between / 
and the depth of the HVF remains even when the 
blobs are crowded as in Model B30. Model B6 still 
has a large space between the blobs, making the 
fraction of small / large. In addition, this model 
is very unlikely to produce the deep absorption 
required to reproduce the objects like SN 2002dj 
and SN 2003du (see. Fig 1 and Fig 4), which needs 

Model B5 (Fig 10, middle) has smaller inter- 
vals between the blobs and a large dispersion in 
the covering factor. This model has a 15% proba- 
bility of showing a deep absorption (/ > 0.9) while 
the fraction of lower / is smaller than in any of the 
models discussed before. The fact that HVFs ex- 
ist in almost all early spectra may suggest that a 
clumpy structure like Model B5 is a possible aver- 
age geometry. 

If we cram more blobs (Model B30; Fig 10, bot- 



tom), the effect of different lines-of-sight is aver- 
aged and the variation disappears. This result is 
similar to that of the ID test with 1 — exp(— r) = 
0.5. Therefore wc can conclude that this geometry 
is not typical of the majority of SNe la, although 
it may exist. 

4.2. Torus Models 

We distinguish a model that has a torus-like 
density enhancement from blob models because 
the origin of the enhancement is expected to be 
interaction with CSM or an accretion disk in this 
case. The main problem with the hypothesis that 
HVFs come from CSM interaction is the fact that 
extremely high mass loss rates are required to 
place the material just outside the explosion. One 
expects that this may be solved if the progenitor 
WD is surrounded by an accretion disk and the 
ejecta collide with the disk, leading to density en- 
hancements on the plane of the disk. Figure 11 
shows the distribution of covering factors of mod- 
els with torus-like enhancements (Model T2, T4 
and T6). These models have a disk-like enhanced 
region with opening angles of 20°, 40° and 60°, re- 
spectively. As shown in the left panels of Figure 
1 1 , the 60° disk has a thickness comparable to the 
diameter of the photosphere. 

In Model T2 the disk is too thin for the cover- 
ing factor ever to become large (top panels in Fig. 
11). If the thickness increases to 40° (Model T4; 
middle panels in Fig. 11), the fraction of "weak" 
HVFs decreases while those of "strong" and "mid- 
dle" HVFs increase. If the the disk is thick enough 
for the photosphere to be totally concealed, the 
fraction of "strong" HVFs increases (Model T6; 
bottom panels in Fig. 11) and the distribution 
of covering factors becomes broad. Only Model 
T6 may produce the observed distribution of HVF 
strength. 

This implies that a very thick disk is required, 
with thickness comparable to the diameter of the 
WD. The material in the disk should be swept by 
the SN ejecta in a radial direction leaving an im- 
print of the disk's angular size after the explosion. 
While this may be a realistic model, the absence 
of hydrogen lines in the early spectra remains the 
strongest argument against the disk origin of the 
HVFs. This problem can be overcome by assum- 
ing low temperature (T ~ 4, 500K) at the shocked 
ejecta (Gerardy et al. 2004), where the hydrogen 
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lines arc less active than metal lines. And see the 
consideration made in Mazzali et al. (2005a) on 
the role of hydrogen in making the Call IR triplet 
broad. They estimate M(H)~ O.OO4M0 in spheri- 
cal symmetry. 

5. Discussion 

The strength of the HVFs {Snv) was parame- 
terized via the line optical depth of the Call IR 
triplet in §2. The observed range of HVFs can 
be reproduced using only geometrical effects if we 
assume that the optically thick region distributes 
discretely. The covering factor that represents how 
the photosphere is concealed by the optically thick 
region acts exactly as 1 — exp(— r) in §2. Among 
the various geometries discussed in the previous 
section, Models Bl and B2 have a very large frac- 
tion of lines-of-sight with / < 0.4, which makes 
almost no HVFs. Since observations suggest that 
a considerable fraction of the earliest spectra show 
HVFs, these models can be ruled out as a standard 
configuration. 

As the interval between the blobs becomes 

smaller, the fraction of lines-of-sight with / < 0.4 
decreases (e.g., from Model B6 to Model B30). 
However, too crowded structures like Model B30 
do not give rise to a variation in the covering factor 
and so the strength of HVF is constant. Model B5 
has a wide range of covering factors from / ~ 0.1 
to / ^ 1. If such a structure is formed from 
the explosion, it can explain the observed variety 
naturally. 

It should be noted that models with a hole can 
be tuned to give the observed distribution of HVF 
more easily than models with "blobs" by changing 
the covering factor, i.e., the size of the hole. How 
such geometries might be produced is however un- 
clear, as they still require a dense high velocity 
region and a line of sight that is never far from 
the hole, because otherwise most SNe would show 
strong HVFs. 

A torus model may reproduce the observed dis- 
tribution of HVF if the disk is geometrically thick 
(~ 60°). Such a thick disk-like enhancement could 
be produced if the accretion disk is thick enough 
to surround the WD before the explosion. The 
opening angle of the accretion disk is likely to be 
reflected on the angular size of a density enhance- 
ment. 



A torus in Model T2 or T4 is so thin that it 
provides a large fraction of "weak" features, which 
are not frequently observed. However, it may be 
possible to reduce the fraction of / < 0.45 if the 
interaction with CSM resulting from spherical or 
bipolar WD wind occurs. Although many addi- 
tional parameters should be introduced to investi- 
gate the possibility of this scenario, this scenario 
seems to work qualitatively. 

When the optically thick region at high veloc- 
ity is introduced as a density enhancement, other 
features are also affected. The Sill A6355 line is 
the most notable line, and its behavior has been 
thoroughly examined (e.g., Benetti et al. 2005). 
Figure 12 shows spectra around Sill A6355. If the 
photosphere is totally covered by dense blobs, the 
absorption minimum becomes u ~ 22, 500km s~^. 
Such high velocity features have not been ob- 
served. Also, the velocity of Sill A6355 is not al- 
ways correlated with that of the Call IR triplet in 
the observed spectra. For example, SN 2002bo has 
a higher Sill A6355 velocity ^ 15, 500km s^-^(at - 
8 days after the maximum) than that of SN SN 
2003kf ~ 12,000km s'^at -9 days) but the Call 
HVF of SN 2002bo {v ~ 22,000km s^^) is slower 
than that of SN 2003kf {v ~ 23,500km s'^). 
Therefore the assumption that HVFs in the Call 
IR triplet are due only to density enhancements 
may not be correct. However there are some SNe 
that have an extended blue wing in Sill A6355. 
This may indicate some degree of density enhance- 
ment also occurs. The combination of density en- 
hancement and abundance enhancement may pro- 
vide a proper correlation. Such a study, however, 
is beyond the scope of this work. It is worth not- 
ing that the mass of additional material is reduced 
from ^ O.lAfQ in the spherical case to ~ 0.05Mq 
in Models B5, B6 and B30. 

The velocity range of HVFs (fnv) is found to 
be the same in all the models and all the lines-of- 
sight as long as the position of the enhancements is 
fixed. If the covering factor is identical, one might 
expect that the position of the absorption mini- 
mum depends on the line-of-sight. However, we 
found that this effect is very slight. The strength 
of the photospheric absorption Sph is not corre- 
lated with the covering factor and it is not affected 
even if the line-of-sight is varied. 
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6. Conclusions 

We tested the properties of HVFs of Call IR 
triplet through ID and 3D modehng. Parame- 
terized ID simulations were used to extract the 
main features that determine the HVF. We defined 
three parameters that govern the HVFs. These are 
the strength of the HVF 5'hv, its velocity range 
Vuv, and the strength of the photospheric com- 
ponent Sph- As for 5'HV) it was shown that ge- 
ometrical effects on models based on high veloc- 
ity blobs or a thick torus can provide enough va- 
riety to cover all the observations. On the con- 
trary, vhv is not affected by line-of-sight effects. 
Therefore we speculate there may be a number of 
blobs that have different velocities. This might be 
possible if the WD rotates and the clumpy struc- 
ture formed by the deflagration flame has different 
properties depending on direction. Such a variety 
is not produced in a torus model. However, it 
may be possible that different disks lead to differ- 
ent enhancements from SN to SN since the degree 
of density enhancement and its velocity range de- 
pend strongly on the density, the total mass and 
radial extent of the disk. 

The strength of the photospheric absorption 
(Sph) may be influenced by the temperature and 
the abundance. Although wc changed this value 
by changing the abundance in this paper, it should 
be studied whether there is any correlation be- 
tween the temperature and Sph- This requires 
knowing the exact temperature and photospheric 
velocity by fitting observed spectra, which will be 
the subject of future work (M. Tanaka et al. in 
preparation) . 

The models that best explain the statistics of 
the observed strength of the HVFs are a blob 
model (B5) with ~ 5 blobs of 80° and a thick 
torus of opening angle ~ 60°. While a torus model 
sounds appealing, the main argument against it is 
the absence of hydrogen lines (but see Gerardy et 
al. 2004, and Mazzali et al. 2005a) and the vari- 
ation of the HVF velocity. A blob model may re- 
sult from the explosion if the mushroom structure 
of the deflagration is not completely washed away 
in a delayed detonation (Gamczo ct al. 2005). 
Although this does not mean that all the SNe 
have a single geometry, the average structure may 
have a distribution of high velocity material sim- 
ilar to that of Model B5. To verify this conclu- 



sion, more early phase observations before maxi- 
mum are needed. A statistical study can constrain 
the structure more accurately. 

Model B5 is different from the 3D structure that 
was obtained for SN 2001el by Kasen et al. (2003). 
Their model has one large blob and an aspherical 
photosphere, and is able to reproduce the observed 
polarization flux. Our results suggest that such a 
geometry cannot be standard (see Model Bl in §4). 
It is completely unknown whether only polarized 
SNe la have such a geometry, consisting of one or 
two blobs in the outermost region. Rotation of the 
WD may lead to suc;li a configuration. Although 
a statistical study of polarization is not available 
yet, it is potentially useful to verify the geometry, 
including the abundance distribution. In addition, 
the statistical study of polarization may be able to 
distinguish the torus models and the blob models 
because the torus-like enhancement tend to pro- 
duce a large polarization (Kasen ct al. 2003). 

Detailed hydrodynamic calculations of 3D de- 
flagration flames suffer from shortage of high ve- 
locity material. It is unknown whether such a re- 
sult indicates our ignorance about the deflagration 
flame or the existence of a transition from a defla- 
gration to a detonation. It is therefore quite im- 
portant to define a suitable explosion mechanism 
from the observational data. This will be achieved 
via a statistical study presented in this paper when 
larger number of early spectra are available. The 
importance of early observations of SNe la cannot 
be understated. 
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(Japan Society for the Promotion of Science) Re- 
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Fig. 1. — Line profiles of the Call IR triplet in the 
early spectra of 4 SNe. Shown are SN 1999ee (9 
days before maximum, Hamuy et al. 2002, red); 
SN 2002dj 11 days before maximum, Kotak et al. 
2005, green); SN 2003cg (8 days before maximum, 
Elias-Rosa et al., in preparation, blue); and SN 
2003du (11 days before maximum, Stanishev et 
al., in preparation, pink). The flux was adjusted to 
have the same continuum flux at A 7500 — 8500A 
by multiplying by a constant. While the HVFs in 
SN 2002dj and SN 2003du are deep, that of SN 
1999ee is weak. SN 2003cg has almost no high 
velocity absorption. 
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Fig. 2. — The effect of a variation in the strength 
of high velocity absorption S'hv on the Call 
IR triplet profile. The line optical depth is 
parametrized as 1 — exp(— r). The various pro- 
files, going from deeper to shallower absorption, 
have 1 - exp(-T) = 0.9,0.7,0.5,0.3 and 0.1, re- 



spectively, at w - 22, 000 - 29, 000km s'^. A vari- 
ation of these parameters does not affect the pho- 
tospheric absorption at all. 
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Fig. 3. — The effect of a variation in the strength 
of photospheric absorption 5ph on the Call IR 
triplet profile. The strength is changed through 
the abundance of Ca. The various profiles, go- 
ing from deeper to shallower, correspond to a 
mass fraction X{Ga) = 0.068,0.034,0.017,0.008 
and 0.004, respectively. The high velocity absorp- 
tion is not affected at all because we fix the value of 
1 — exp(— t) = 0.7 and the velocity separation be- 
tween the photospheric component and the HVF 
is too large. 
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Fig. 4. — The observed Call IR triplet (solid lines) compared with synthetic spectra computed by adjusting 
all three parameters (S'hv,whv and Sph- AH spectra are normalized to their continuum level. The red 
dashed lines are the line profiles computed including the high velocity components, while the blue dotted 
lines are the profiles computed with only the photospheric component. S'hv corresponds to 1 — exp(— r), 
where r is the line optical depth of the Call IR triplet in the high velocity region. Spi^ is varied through 
the abundance of Ca [X(Ca)]. SN 1999ee (9 days before maximum, -9days; top left) is reproduced with 
1 - exp(-r) = 0.5, whv ^ 22,000 - 27,000km s"\ and X{Ca) = 0.004. SN 2002dj (-11 days; top right) is 
reproduced with l-exp(-T) = 0.7, urv ~ 22, 000 - 29, 000km s-\ andX(Ca) = 0.008. SN 2003cg (-8 days; 
bottom left) is reproduced with 1 - cxp(-T) = 0.2, vhv ~ 22, 000 - 25, 000km s-\ and X{Ca) = 0.004. SN 
2003du (-11 days; bottom right) is reproduced with 1 - exp(-T) = 0.9, whv ~ 22, 000 - 25, 000km s"\ and 
X{Ca) = 0.008. 
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Fig. 5. — The geometry of Model Bl (left) and synthetic spectra obtained viewing this geometry ft-om 
different orientations (right). The opening angle of the blob is ~ 80°. The red dashed lines in the right 
panel show the synthetic spectra corresponding to viewing angles of 0°, 21°, 33° and 42°, respectively, going 
from deep to shallow absorption. The angle is measured from z-axis (the direction of the blob). For angles 
closer to the equatorial plane the depth of the HVF is almost the same as in the spherical model without 
any enhancement (thick line). 
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Fig. 6. — 3D synthetic spectra computed from 
Model Bl (red) compared with the ID results of 
Figure 2 (blue). While the ID results have been 
obtained for shells with different opacities, the 3D 
spectra correspond to the same model viewed from 
different angles, and therefore differ by the cover- 
ing factor /. The black line is the model without 
high velocity enhancement, or / = in 3D. The 
ID spectra have 1 — exp(— r) ~ 0.9, 0.7, 0.5 and 
0.3, from deeper to shallower. The 3D spectra 
have covering factor / = 1.0,0.88,0.66 and 0.46, 
from deeper to shallower. ID and 3D results are 
related. The difference between the values of / 
and 1 — exp(— r) comes from the fact that the blob 
in the 3D model is not infinitely thick (see §4). 
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Fig. 7. — Probability distribution of the covering 
factor / in Model Bl. The horizontal axis shows 
the covering factor and the vertical axis shows 
the frequency of the covering factor if we see this 
model from various lines-of-sight. 
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Fig. 8. — Probability distribution of the observed 
HVFs. 8 SNe from Mazzah et al. (2005 and 2005b) 
are distributed to 4 bins by fitting the depth of the 
HVFs. Note the number of the bin is not identical 
with the other histogram of the covering factor /. 
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eometry of Model B2 (left) and the distribution of the covering factor (right). 
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Fig. 10. — Top: the geometry of Model B6 (6 blobs with opening angle 60deg) and its distribution of covering 
factor. Middle: same as Top but for Model B5 (5 blobs with opening angle 80deg). Bottom: same as Top 
but for Model B30 (30 blobs with opening angle SOdeg) 
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Fig. 11. — Same as Figure 9 but for torus Model T2 (opening angle 20 deg, top), Model T4 (opening angle 
40 deg, middle) and Model T6 (opening angle 60 deg, bottom). 
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Fig. 12. — Synthetic spectra around Sill A6355. 
The thick Una is a ID spectrum without any en- 
hancement. Dashed hnes are 3D spectra with 
covering factor / = 1.0,0.88,0.66 and 0.46 from 
deeper to shaUower. If the photosphere is totally 
hidden (/ ~ 1), the feature forms at a high veloc- 
ity that has never been observed in any SN. 
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